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We present a calculation of the antineutrino flux produced by the reactors at the Laguna Verde Nu¬ 
clear Power Plant in Mexico, based on the antineutrino spectra produced in the decay chains of the 
fission fragments of the main isotopes in the reactor core, and their fission rates, that have been cal¬ 
culated using the DRAGON simulation code. We also present an estimate of the number of expected 
events in a detector made of plastic scintillator with a mass of 1 ton, at 100 m from the reactor cores. 
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I. INTRODUCTION 

In recent years, there has been significant progress in 
understanding the properties of neutrinos, most notably, 
since 2012 the results of the reactor neutrino experiments 
Daya Bay [T], Reno and Double Chooz [2] have deter¬ 
mined that the mixing angle Ois is nonzero with a high 
level of significance. The current degree of development 
of neutrino detection technology is close to making a re¬ 
ality to non-intrusively monitor the operational status, 
power level and fissile content of a nuclear reactor in real 
time using detectors placed at distances of a few tens of 
meters. 

Monitoring nuclear reactors through its antineutrino 
flux is a complimentary and promising new tool for super¬ 
vising the operations of nuclear plants, which are bound 
to operate according to protocols established by the In¬ 
ternational Atomic Energy Agency (IAEA). This agency 
is responsible for ensuring that nuclear reactors world 
wide operate legally, preventing the diversion of fissile 
material to activities which could lead to the manufac¬ 
ture of weapons. 

Many research groups worldwide (i-ls!) have studied 
this application aiming at adding antineutrino detection 
to the techniques used to implement reactor safeguards. 
In a medium term plan (5 to 8 years) IAEA [9] has pro¬ 
posed that antineutrino detectors for safeguard applica¬ 
tions should be able to provide information on the ther¬ 
mal power, fissile content and operational status of reac¬ 
tors while deployed aboveground in a compact volume to 
minimize its intrusiveness with the nuclear power plant 
operations. The SONGS I experiment [3] using a Gd 
loaded liquid scintillator detector, deployed at 25 m from 
the reactor core of the San Onofre Nuclear Generation 
Station, successfully proved that the reactor ON/OEE 
cycle and the thermal power stability of a reactor can be 
measured using a compact and simple detector. PANDA 
[4], lead by Tokyo University, used a segmented plas¬ 
tic scintillator detector with modules wrapped with Gd 
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coated sheets in between. This type of detector makes it 
possible to use information of the event topology to tag 
antineutrino events and to discriminate them from back¬ 
grounds such as fast neutrons. Eurthermore, it can be 
transported and operated easily inside a compact vehicle, 
and has the advantage that unlike the liquid scintillators, 
is non-flammable. 


Mexico has two nuclear fission reactors at the Laguna 
Verde Nuclear Power Plant (LVNPP), which operates 
commercially since 1990 (Unit I) and 1995 (Unit 2). This 
is the only nuclear power plant in the country and gen¬ 
erates about 5% of its total electric power production. 
In this work we calculate the antineutrino flux produced 
by these reactors and give some estimates of the number 
of events expected in a generic plastic scintillator detec¬ 
tor located at a distance of 100 m from the reactor cores. 
Knowledge of the antineutrino flux produced by a partic¬ 
ular reactor is an important step towards its monitoring 
using this method. 


While abundant literature exists about the calcula¬ 
tion of neutrino fluxes from pressurized water reactors 
(PWR) [To], much fewer information is available about 
neutrino production at boiling water reactors (BWR) 
(See for example HI]). These two reactor types oper¬ 
ate under the same basic physical principles regarding 
nuclear fission, but differences in their operational pro¬ 
cedures may effect differently the neutrino flux they pro¬ 
duce, as well as its evolution along the reactor fuel cycle. 


This paper is organized as follows: in section [Hj main 
technical data on the LVNPP and its reactors are de¬ 
tailed. In section [lll| a description about the production 
of antineutrinos from nuclear reactors is given, and the 
calculation of antineutrino flux is presented. In Section 


IV parameters and characteristics of the DRAGON sim¬ 


ulation code used to calculate fission rates of four major 
fissile isotopes for this paper are described. Section \V\ 
shows the event rates or interactions that are expected 
to be observed with a I ton plastic scintillator detector 
placed at 100 meters from the LVNPP reactors. Einally, 
sectior fVT] presents the conclusions. 





II. THE LAGUNA VERDE NUCLEAR POWER 
PLANT (LVNPP) 
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The LVNPP is located in the municipality of Alto 
Lucero, Veracruz, Mexico. This plant has two twin units, 
each equipped with a BWR-5 second generation reactor 
with a Mark II containment design, supplied by General 
Electric with a capacity of 2027 MWth, and net electri¬ 
cal output of 805 MWe. Both reactors operate with en¬ 
riched uranium as fuel, and demineralized water as mod¬ 
erator and coolant. It is the only nuclear power plant in 
the country and generates about 5% of the total electric 
power production in Mexico [12]. 

Each nuclear reactor core can be approximately repre¬ 
sented by a cylinder ^4 m in height and 4 m in diameter 
containing 444 fuel assemblies and 109 control bars ar¬ 
ranged as shown in Eig. A fuel assembly is a square 
prism with sides of ^13 cm and 4 m in height, contain¬ 
ing an array of fuel rods (configurations of 8x8, 10x10 
or 12x12 have been either used or considered). Within a 
given fuel assembly some of the rods contain uranium ox¬ 
ide (UO 2 ) with various levels of enrichment, others con¬ 
tain a mixture of UO 2 with gadolinium oxide (Gd 2 03 ) 
at varying concentrations, and some of them are hollow 
in order to allow the flow of the refrigerant through them. 
Combining rods with different levels of enrichment and 
Gd 203 concentrations permits to adjust the properties of 
the medium for the production, absorption and diffusion 
of neutrons differently at various location across the core. 



9 68 natural uranium assemblies (0.711% 

96 half enrichment assemblies (1.76% 

280 high enrichment fuel assemblies (2.19% 

109 control bars 

FIG. 1. Schematic of a cross section of one of the LVNPP 
reactor cores [19]. The ”red” assemblies (0.711%) are located 
on the perimeter of the reactor core, while the ’’blue” assem¬ 
blies (1.76%) are arranged in the central area. There are 109 
control bars represented as crosses in the figure. 


III. DRAGON SIMULATION OF THE 
REACTOR CORE 

We used the numerical simulation code DRAGON m 
with the modifications described in m to implement a 
simulation for the LVNPP reactors and extract the time 
evolution of the fission rates (number of fissions per sec¬ 
ond) of each of the main fisile isotopes in the core. Several 
research groups worldwide have implemented this simula¬ 
tion for PWR reactors, aiming at the same objective m- 
The code solves the neutron transport equations of in¬ 
dividual fuel assemblies, whose detailed composition is 
given as an input. The overall behaviour of neutrons 
across the core can be calculated by summing over the 
various types of assemblies that compose a given core 
configuration. 

As described in section |TI| a geometry with 444 fuel 
assemblies with 3 different types of enrichment was 
considered. The initial fuel load that was assumed in this 
work is that shown in Fig. Note that in this model 
no assembly contains plutonium at the beginning of the 
operation cycle. 

In Fig. we illustrate a single fuel assembly (in this 
case with a 10x10 rod configuration) corresponding to 
one with average enrichment of 1.76% (blue assembly). 
The dimensions of a single unit cell, (the place occupied 
by a single fuel rod) are shown in Fig. The three 


types of assembly (red, yellow and blue) share the same 
geometry, however, the composition the fuel rods within 
a given assembly is, in general, different. The number of 
assemblies of each type, Na is: Nred = 68, Ny^uow = 280 
and Nuue = 96, totaling 444. 

We considered eight different fuel rod compositions, la¬ 
beled CX, with X = 0,1, 2, 3,4, 5, 6, G, and used them to 
construct the various assemblies with the desired en¬ 
richment and Gadolinium concentration. Only fuel rods 
of type G contain Gadolinium. Table |T] shows the number 
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FIG. 2. Arrangement of fuel rods for blue assembly (half 
enrichment). Left: DRAGON code output. Right: Bar type 
(see Table [lT|. Notice that initially there are no fuel bars with 
Pu. 
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FIG. 3. Cross section of a single unit cell. A fuel bar (in red) 
with a zircaloy cladding (white) immersed in demineralized 
water as moderator. A complete assembly is formed by 92 
of these cells. The values for Rl, R2 and L are described in 
Table [ml 

of fuel rods of each type that are used to build each of 
the three assemblies. As an example, Table [TT| shows how 
we constructed the blue assembly. 


Assembly 

GO 

G1 

G2 

G3 

G4 

G5 

G6 

GG 

Red 

92 

0 

0 

0 

0 

0 

0 

0 

Blue 

59 

0 

0 

3 

4 

6 

8 

12 

Yellow 

34 

9 

19 

0 

0 

0 

0 

30 


TABLE I. Number of fuel bars of type CX for each assembly. 
There is a total of 92 bars for each assembly. 

Table m lists the values used for the various parame¬ 
ters required by the simulation. Tests showed a small de¬ 
pendence of the fission rates on variations of the coolant, 
moderator and fuel temperatures, while maintaining the 
thermal power constant. 

A simulation was run for each of the three types of 
assembly (red, blue or yellow in Fig. [^, providing the 
fission rates for each type of fuel rod, and for 

each one of the fissile isotopes i = 1,2, 3,4 (^^^U, 

^^^Pu and ^^^Pu, respectively). 

To obtain the fission rates for a given assembly, //^, we 
multiply the fission rate for each type of fuel rod times 
the number of rods of its type in the assembly, N^ods^ 
and sum over all fuel rod types in the assembly, i.e: 

/^= E Nrods{CX-A)ff^''^^ (1) 

cxeA 

Finally, to calculate the fission rates of isotope i in the 
complete reactor core, fi (Fig. |^, we multiply the num¬ 
ber of assemblies of each type A, times the corresponding 
fission rates for that assembly: 

A 


Rod type 

# of rods 

2-64^ 

(%) 

235 

(%) 

238 -^ 

(%) 

016 

(%) 

Gd 

(%) 

GO 

59 

0.005 

0.627 

87.517 

11.852 

- 

G3 

3 

0.029 

2.468 

85.651 

11.852 

- 

G4 

4 

0.029 

2.821 

85.298 

11.852 

- 

G5 

6 

0.029 

3.173 

84.945 

11.852 

- 

G6 

8 

0.029 

3.482 

84.637 

11.852 

- 

GG 

12 

0.028 

3.142 

79.691 

11.933 

5.206 


TABLE 11. Isotopic composition of fuel bars for blue assembly 
in Figg UO 2 and U 02 -Gd 203 bars with an overall average 
enrichment of 1.76% of (relative to the total U). Enel rods 
with Gd admixture have an enrichment of 5.2% (in total) and 
3.792% (relative to U). 


The simulation was run simulating time intervals of 5 
days of evolution of the nuclear reactor core until com¬ 
plete a full operation cycle of 400 days. 

As is expected for an initial load without plutonium, 
the fission rates of ^^^Pu and ^^^Pu start at zero and 
rapidly grow to become comparable to those of the U 
isotopes within a few weeks of operation. A more realistic 
situation should consider that between two fuel reloads, 
there will always be some remnant of plutonium inside 
the reactor, because when this action is performed, only 
a fraction of the spent fuel is replaced by new material. 


Thermal Power (of each reactor) 

2027 MWth 

Specihc power (of each reactor) 

20.43 MW/Ton 

Moderator 

Demineralized water 

Temperature 

600 K 

Density 

0.720 g/cm® 

Goolant 

Demineralized water 

Temperature 

400 K 

Density 

0.720 g/cm® 

Enel assemblies 

444 

Red assemblies 

68 

Blue assemblies 

96 

Yellow assemblies 

280 

Single unit cell (LxL) 

1.295 cmx 1.295 cm 

Enel bars per assembly 

92 

Gomposition 

U02/U02-Gd203 

Gladding 

Zircaloy 2 

Zircaloy density 

5.821 g/cm^ 

Zircaloy temperature 

600 K 

Enel temperature 

900 K 

Enel density 

10.079 g/cm^ 

Enel bar radius (Rl) 

0.438 cm 

Zircaloy cladding radius (R2) 

0.513 cm 

Enel bar height 

400 cm 

Goolant pipes per assembly 

2 

Gomposition 

Zircaloy 2 

Pipe radius (occupied by water) 

1.0 cm 

Pipe radius (occupied by zircaloy) 

1.2 cm 


TABLE III. Parameters of the BWR core used in the simula¬ 
tion with DRAGON code. 
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Days 


FIG. 4. Fission rates (normalized) obtained with DRAGON 
for the 4 main fissile isotopes in the 444 fuel assemblies that 
make up the core of a BWR-5 of 2.027 GWth like the LVNPP 
reactor. 


IV. THE ANTINEUTRINO FLUX 

The power output of a nuclear reactor originates from 
the energy released in the fission of heavy elements, such 
as uranium and plutonium, into lighter fission fragments 
which are often unstable. The beta decay, (A, Z) 

(A, Z + 1) + e“ + Pe, of the fission fragments produces a 
large number of electron antineutrinos Pe which are emit¬ 
ted isotropically from the core. During the fuel ’’burn¬ 
ing” process, uranium isotopes breed ^^^Pu and ^^^Pu. 
The latter may in addition to the undergo a fission 
process only with slow (thermal) neutrons, while is 
fissionable by fast neutrons only. The decay of ^^^Pu pro¬ 
duces substantially less antineutrinos than the decay of 
235 u in the same energy range, therefore during a typical 
reactor fuel cycle, the amount of antineutrinos decreases 
as uranium content decreases and the concentration of 
plutonium increases. 

Each fission releases on average ^200 MeV of en¬ 
ergy and produces ^6 Pe 3 beta decays per fission 
fragment), with energies below ^10 MeV. This sets the 
number of antineutrinos emitted by a typical reactor to 
^ 2 X 10^^Pe/«5 per GWth of thermal power. The energy 
spectrum of the antineutrinos depends on the fuel com¬ 
position at a given time (^^^Pu antineutrinos are slightly 
less energetic than those of fission products). 

The production of reactor antineutrinos is not exclu¬ 
sively through beta decay of the fission fragments of 
the four main fissile isotopes (^^^U, ^^^Pu and 

^^^Pu).The neutron capture in ^^^U(n, 7 ) 
crates this particles in a process that contributes about 
17% of the total antineutrino flux. This process occurs 
when a nucleus of captures a neutron, leading to 
the following reaction: Np ^^^^Pu, 

and produces two antineutrinos through two beta decays 
of This process occurs at a rate of ^0.6 per fission, 
and produces antineutrinos with energies below 1.3 MeV. 

There is no established experimental method for mea¬ 


suring the flux of very low energy antineutrinos. For 
a technological application such as reactor monitoring, 
only antineutrinos with energies above the threshold of 
the Inverse Beta Decay (IBD) reaction on free protons 
(Ejy >1.8 MeV) can be realistically considered. IBD in¬ 
teractions of neutrinos on free protons can be identified 
through the correlation of the prompt positron signal and 
a delayed neutron capture signal. In general, the spec¬ 
trum of the antineutrinos from a nuclear reactor detected 
via IBD has a peak between 3.5 and 4 MeV. Since the 
number of emitted antineutrinos and their average en¬ 
ergy depend on the amounts of ^^^Pu and 

^^^Pu isotopes present in the core, the measured energy 
spectrum provides a direct image of the fuel composition 
of the reactor core at a given time. 

Let us now consider the antineutrinos emitted by one 
of the reactors at the LVNPP. Ignoring the neutron cap¬ 
ture contribution, the antineutrino flux above 1.8 MeV 
can be calculated from the fission rates /i(t), and the 
antineutrino energy spectrum Si{Ep) of each isotope i, 
where i = 1, 2, 3,4 corresponds to each of the main fissile 
isotopes of the reactor core: ^^^Pu and ^^^Pu. 

Assuming neutrinos are emitted isotropically, using a 
sphere of radius i7, the flux at a time t for each isotope i 
is calculated as: 

(3) 

with R the distance from the reactor core to the detector. 

In our calculation of the antineutrino flux, we have 
used the new predictions of the energy spectra includ¬ 
ing the effect of the reactor antineutrino anomaly, as re¬ 
ported by m, m and m- The antineutrino energy 
spectra Si{Ey) between 1.8 and 8 MeV are shown in Eig. 

Eor the purposes of reactor monitoring, the effect of 
this anomaly can be ignored as long as only comparisons 
between the relative changes of the spectrum throughout 
the reactor operation cycle are made. 



FIG. 5. Energy spectra of the emitted antineutrinos per fis¬ 
sion of ^^®Pu and ^^^Pu. The dotted line at 1.8 

MeV indicates the threshold for the Inverse Beta Decay (IBD) 
process. 
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FIG. 6. Antineutrino flux produced by the fission of the 4 
main fissile isotopes at 100 meters from one of the LVNPP 
reactors of 2.027 GWth at day 5 (up) and 200 (bottom) of 
operation. The gray curve (on top) is the total flux. 

Fig. shows the antineutrino flux for each of the main 
fissile isotopes of one of the reactors in the LVNPP, at day 
5 (top) and then at day 200 (bottom), in the latter the 
flux can be considered as stable. The total flux (the sum 
of the four isotopes) for each time period is also shown 
in the gray curve in this figure. The total flux at day 200 
is shown in table lEl 

Two systematic uncertainties are associated with the 
calculation of the antineutrino flux: the error due to the 
energy spectra of antineutrinos S{S{Ejy)), and an error 
6{Pth) associated with fluctuations in the thermal power 
of the reactor core, assumed to be 5% along one fuel cycle. 
The total error is the sum in quadrature of these uncer¬ 
tainties. The uncertainties in the antineutrino spectra 
are taken from Refs. El and m, and are of the or¬ 
der of 2-5% for the dominant parts of the spectra from 
235u 239p^^ and ^^^Pu, and between 10-17% for the dom¬ 
inant part of the spectrum from 

The total flux per reactor at day 200 from our calcula¬ 
tion is compared in Fig. [^to that obtained according to 
the parametrization suggested by Vogel and Engel [22] . 
setting the fission rates of the four main fissile isotopes to 
these output by the simulation (2.7 xlO^^ fis/s for 
7.7 xlO^^ fis/s for 2.8 xlO^^ fis/s for ^^^Pu, and 

1.5 xlO^^ fis/s for 241 Pu). The two calculations agree 


E [MeV] 

^totai [MeV-^ s-%m-^] 

(5(sist) [%] 

1.800 

6.921E+10 

7.2 

2.048 

6.192E+10 

3.6 

2.296 

5.185E+10 

3.6 

2.544 

4.304E+10 

3.7 

2.792 

3.634E+10 

3.7 

3.040 

3.052E+10 

3.7 

3.288 

2.542E+10 

3.7 

3.536 

2.066E+10 

3.8 

3.784 

1.671E+10 

3.8 

4.032 

1.336E+10 

3.9 

4.280 

1.045E+10 

4.0 

4.528 

8.026E+09 

4.0 

4.776 

6.147E+09 

4.1 

5.024 

4.841E+09 

4.2 

5.272 

3.761E+09 

4.3 

5.520 

2.871E+09 

5.2 

5.768 

2.157E+09 

5.4 

6.016 

1.584E+09 

5.5 

6.264 

1.179E+09 

5.6 

6.512 

8.900E+08 

5.8 

6.760 

6.692E+08 

5.9 

7.008 

4.752E+08 

7.0 

7.256 

3.027E+08 

7.4 

7.504 

1.881E+08 

7.8 

7.752 

1.210E+08 

8.4 

8.000 

6.763E+07 

9.8 


TABLE IV. Total antineutrino flux at day 200 of one of the 
reactors of the LVNPP, at a location of 100 m from the reactor 
core. 


well within the assumed uncertainties. 

Spent fuel is stored temporarily in a nearby place to 
the reactor core and has a very small contribution to the 
antineutrino flux below the IBD reaction threshold, and 
was not considered for the calculation in this work. 


V. EVENT RATES 

Although the IBD reaction has a very small cross sec¬ 
tion (^10“^^cm^), the enormous flux emitted by a nu¬ 
clear reactor allows the signal to be observed with a rela¬ 
tively small detector located a few tens of meters from the 
reactor core. The total IBD cross-section as a function 
of the neutrino energy is shown in Eig. (right scale). 
Here we consider 1 ton of polyvinyl-toluene (PVT) plas¬ 
tic scintillator located at 100 m from each of the reactor 
cores, similar to the detector considered in [4]. 

The number of IBD interactions associated with neu¬ 
trinos from isotope i, emitted with energies between E^in 
and Ejnax^ and occurring in a time interval At (AVev(i)) 
is calculated by integrating the product of the flux the 
IBD cross section (j{Ep) [23], and the number of targets 
(in this case free protons) Vp, over the energy interval: 
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FIG. 7. Total antineutrino flux per reactor at day 200 (solid 
curve with error band). The error includes the sum of the 
relative uncertainties of each of the 4 main fissile isotopes in 
the reactor. The dashed curve shows the parametrization sug¬ 
gested in with the fission rates calculated by DRAGON 
for the LVNPP. The continuous line shows the IBD cross sec¬ 
tion as given by [23] 



FIG. 8. Number of events accumulated over 2 days of expo¬ 
sure along a complete cycle of 400 days at LVNPP with two 
reactors in 1 ton of plastic scintillator. Gurves for detection 
efficiencies (lower to higher) of 5%, 20%, 30%, and 100% are 
shown 


A V /-x pEmax 

= Np (4) 

Depending on the characteristics of the detector, only 
a fraction of these interactions will be observed. For a 
simple estimate, let this fraction be a uniform efficiency 
factor e. For neutrinos coming from nuclear reactors, 
Emin is the threshold of 1.8 MeV and Emax is ^ 10 MeV. 
The number of detected events for each isotope is then 
approximated by: 


of ^365 events/day, operating a 360 kg detector with 
the reactor OFF for 21 days. Scaling to a 1 ton detector 
corresponds to ^1013 events/day. Since our expected 
rate at LVNPP is 26 events/day at the location of 100 m 
from the reactor cores, a background measurement at the 
level of 0.5% (±5 events/day) would be necessary for a 5cr 
detection of the antineutrinos, comparing the ON/OFF 
event rates with a similar detector. 


VI. CONCLUSIONS 


A/ P-^max 

7Ve,(i) = e X SiiE,)Mt)a{E,) dE„ (5) 

where At is a small exposure time compared with the 
time scale in which fission rates fi{t) change appreciably. 
Calculation for longer times requires an integration over 
time. 

Finally the total number of events is equal to the sum 
of the events for each isotope: 

4 

(total) ^ ^ ^ev{i) (0) 

i=l 

These are shown in Fig. considering three differ¬ 
ent detector efficiencies (5%, 20%, 30% and 100%) and 
counting for two consecutive days. 

The authors of |4| estimate an antineutrino detection 
efficiency of 9.24% and an antineutrino event rate of 147 
events/day in a 1 ton detector with an uncertainty of the 
order of 30%. This is consistent, within errors, with our 
expected rate of 280 events/day (stable after day 350), 
once scaling for reactor power, distance, and assuming 
the same efficiency. The authors measured a background 


We calculated the antineutrino flux produced by the re¬ 
actors at the LVNPP in Mexico using the DRAGON sim¬ 
ulation code in the energy range 1.8 MeV<Ejy<S MeV. 
The simulation provided the time evolution of the fis¬ 
sion rates and hence the change of the flux along the 400 
day long fuel cycle, starting with an initial load with no 
plutonium. A total systematic uncertainty on the flux 
ranging from 3.6% to 9.8%, depending on the antineu¬ 
trino energy, was calculated assuming a 5% uncertainty 
on the reactor power, and the uncertainties reported for 
the antineutrino spectra of the dominant fissile isotopes 
( MM) .Our calculation agrees well with a frequently 
parametrization [22| when the fission rates output by the 
simulation of the LVNPP BWR reactors are used. 

The simulation code required the knowledge of generic 
parameters of the BWR reactors, at LVNPP, available 
in the literature, and showed the fission rates to have a 
moderate dependence on variations of the coolant, mod¬ 
erator and fuel temperatures for a fixed thermal power. 

A segmented plastic scintillator detector of the type 
considered in |4] located at 100 m from the two reac¬ 
tors could detect the neutrinos emitted, provided a back¬ 
ground measurement with the reactors OFF is performed 
at the 0.5% level. 
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